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P/ Abstract 



We show that "top-down" mechanisms of Ultra-High Energy Cosmic Rays 



■ which involve heavy relic particle-like objects predict Galactic anisotropy of 

> 

highest energy cosmic rays at the level of minimum ~ 20%. This anisotropy 

OO 

is large enough to be either observed or ruled out in the next generation of 

o ' 

00 ■ experiments. 
On. 

The origin and nature of Ultra-High Energy Cosmic Rays (UHE CR) with energies 

a: 

above 10 eV [|I| is one of the actively debated issues in modern astrophysics. At such a 
high energy the mean free path of protons, nuclei and photons is much shorter than the size 



of the Universe due to their interaction with cosmic microwave and radio background 



Protons with E > Eqzk ~ 5 x 10 19 eV can only come from distances less than Rqzk ~ 50 
Mpc; the corresponding region for nuclei and photons is even smaller. The energy spectrum 
of cosmic rays is thus expected to have a rapid falloff at energies E ~ Eqzk (the so-called 
GZK cutoff). 

The observation of few events with energies exceeding 10 20 eV seems to indicate that 
noticeable fraction of UHE CR comes from relatively nearby (R < Rqzk) sources. The latter 
is not easy to reconcile with conventional astrophysical mechanisms Q as at these energies 
and distances the cosmic rays are not substantially deflected by magnetic fields and point 
in the direction of their source ||. Corresponding astrophysical sources have not yet been 
identified. This prompts to consider particle physics mechanisms of UHE CR production, 



which are usually referred to as "top-down" mechanisms. Several such mechanisms involving 
topological defects and decays of primordial heavy particles 0-0] have been proposed. 
In this Letter we show that large class of the "top-down" mechanisms predict significant 
excess of highest energy events in the direction towards the center of our Galaxy. Thus, 
these mechanisms can be tested experimentally by studying the asymmetry in the angular 
distribution of UHE CR. The predicted asymmetry is large enough to be either observed or 
excluded by future experiments |Tl| . 

Our consideration concerns mechanisms in which production of UHE CR involves heavy 
relic particle-like objects which behave as part of cold dark matter (CDM). We call such 
mechanisms CDM-related. The crucial feature of these mechanisms is that the distribution of 
sources of UHE CR in the Universe between galaxies and intergalactic space is proportional 
to that of CDM (cf. ref. fl2f ) and does not depend on their nature. In particular, the sources 
are mainly concentrated in galactic halos, so that their average densities in the Universe (n) 
and in galactic halo (n h ) are related by 

n ^CDMPcrit ln _5 /1N 

— ~ — = ~ 10 . (1) 

nh Phalo 

On the contrary, the distribution n(x) of the sources in the galactic halo depends on their 
interaction with each other and with other matter and does not necessarily follow that of 
CDMQ. For the decay-type mechanisms n(x) is simply the density of decaying particles. 
For mechanisms based on collisions (see e.g. JT3]) n(x) oc n 2 (x), where n(x) is the density 
of colliding particles. Note that in the latter case the distribution n(x) is typically more 
concentrated around the galactic center. 

The observed flux of UHE CR can be divided into Galactic and extragalactic parts, 

J Jcxt ")~ J hi 

where 
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h = (' I ^n(x) (2) 



halo X 2 

is the contribution of our Galaxy and 

j ext = C 4nR ext n 

has extragalactic origin. Here R ex t — -^Universe ~ 4 Gpc for energies below Egzk and. 
-Rext ~ 50 Mpc for energies above -Egzk- The constant C is the same in both equations. 
Following ref. || we note that eq.(|l]) allows to estimate the relative magnitude of these two 
contributions, 

Jext -Rext 71 -Rext 1 n— 5 /o\ 

— = o>——— ~ a —fr~ x 10 > (y) 

jh Rh n h R h 

where Rh ~ 100 kpc is the size of the Galactic halo and a is the constant of purely geometrical 
origin, 



3 Ir<R h d 3 



xn\x) 



a = " Jr<Rh ^-v-^ _ (4) 

R/i Ir<R h ~Jf n (. X ) 

Here r(x, 0) = (x 2 + r 2 — 2xr cos 0) 1 / 2 is the distance between current point and the Galactic 
center while r = 8.5 kpc is the distance to the Galactic center. 

Numerical value of a is close to 1 for uniform distribution n(x) = const and can be small 
for distributions concentrated around the galactic center. Although n(x) does not have to 
coincide with CDM distribution in the halo, we consider as two examples the isothermal 



halo model 14 



(r£ + r z ) 



and more realistic distribution of ref. [15 



n(r) oc , (6) 

yJ(r* + r>)(R h + r) 2 

which we have arbitrarily regularized at r = by introducing the core size r c . The value 
of a is a ~ 0.15 and a ~ 0.5 for distributions and (|^), respectively, with no strong 
dependence on r c in the range r c = 2 — 10 kpc. 



From eq.(|3|) we find 

~ a for E < E G zk, 

Jh 

J — ~ 10~ 2 xa for E > E GZK . (7) 

Jh 

Therefore, at E < Eqzk the Galactic and extragalactic contributions can be comparable (al- 
though the Galactic one is probably somewhat larger), while at E > Eqzk the extragalactic 
part is suppressed by a factor ~ 1CT 2 . In either case a substantial fraction of the observed 
UHE CR should come from the halo of our Galaxy. In this respect our conclusions agree 
with that of ref. 0. 

The Galactic part of the total UHE CR flux, jh, is anisotropic due to our position at 8.5 
kpc from the center of the Galaxy. The anisotropy can be obtained from eq.(^), 

jh{6) oc J dxn(r(x,9)). 

Fig. la shows the anisotropy jh(Q) / jh(n) as a function of the core radius for the trial dis- 
tributions fl5|) and (P). For n(x) = const the anisotropy is minimum and constitutes about 
20%. Fig. lb shows corresponding angular dependencies of jh(@) at r c = 5 kpc. As can be 
seen from the picture, the anisotropy of the galactic contribution is at least ~ 20% and can 
be much larger if n(x) is concentrated around the galactic center. Also, it should be noted 
that the anisotropy depends exclusively on n(x) and does not depend on energy since cosmic 
rays with energy E ~ E GZ k are deflected by the Galactic magnetic field by only ~ 3° 

In turn, the extragalactic contribution j ext consists of the isotropic part coming from 
distances R 3> 50 Mpc, which is comparable in magnitude to jh and is present only at 
E < Eqzk-i and the contribution from our "neighborhood" R % 50 Mpc. The latter should 
have peaks in the direction of nearby galaxies and clusters. The contribution of such a peak, 
Sj ext , equals 

6j ext R\ M 

= ot 

Jh 3i? 2 M G ' 

where R is the distance to the astronomical object, M is its mass, and Mq is the mass of 
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our Galaxy including halo. For instance, contributions from Andromeda Nebula and Virgo 
Cluster are comparable and close to 10~ 2 x a, in agreement with eq.(|7|) and ref. ||. 

Since at energies above the GZK cutoff the extragalactic contribution is negligible, non- 
observation of the anisotropy at the level of ~ 20% would rule out the CDM-related mech- 
anisms of UHE CR. The observation of the Galactic anisotropy would allow to reconstruct 
the density profile n(x) and, possibly, the distribution of CDM in the Galactic halo. 

At energies below the GZK cutoff, the anisotropy is smaller due to the relative enhance- 
ment of the isotropic extragalactic part. Since anisotropy does not depend on energy and 
can be measured at E > Eqzk, it is possible, in principle, to determine the magnitude 
of the extragalactic contribution. Provided the CDM-related mechanisms are dominant at 
E < Eqzk and the coefficient a is known, the ratio jh/jext could give, in view of eqs.(|l|) and 
(0), an important information about the distribution of matter in the Universe. 

Current data is not enough to draw definite conclusions about the angular distribution 
of highest energy cosmic rays both because of very limited statistics and the absence of 
data in the South hemisphere where the Galactic center is situated. However, since the 
anisotropy predicted by the CDM-related mechanisms is large, it will be either observed or 
excluded already in the next generation of experiments [IT]. Among these the Pierre Auger 
Project has the best potential due to large number of expected events (600-1000 events with 
E > 10 20 eV in 10 years) and the ability to see both hemispheres. 
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FIGURES 




r c , kpc 9 

FIG. 1. a) The anisotropy jh(0)/jh(^) as a function of the core size r c for the density profiles 
(solid line) and (||) (dashed line), b) The corresponding angular distributions at r c = 5 kpc. The 
dotted line shows the angular distribution for n{x) = const (i.e., when the anisotropy is minimum). 
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